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ABSTRACT 

Aims. Searching for planets in the atmosphere of AGB stars is difficult, due to confusion with the stellar wind and pulsations. The aim 
here is to provide a complementary strategy for planet searches in such a dense environment. 

Methods. The polarization properties of SiO masers, especially their circular polarization, are, under certain conditions, good tracers 
of rapid magnetospheric events. A Jovian planet with a magnetosphere whose dipole axis is misaligned with its rotation axis naturally 
provides such conditions. Here I present several models showing that the polarization will be periodically modulated. 
Results. The linear and circular polarization of an SiO maser in a planetary magnetosphere is modulated by the precessing dipole 
component of the latter. The effect is measurable in saturated masers, while unsaturated masers only exhibit weak changes, because of 
dilution effects, and because the circular polarization there stems from the Zeeman effect making it as weak as for thermal radiation. 
The situation would change if anisotropic pump- and loss-rates were included, which would increase the fractional linear and, via 
magnetorotation, the circular polarization of the modulation. 

Conclusions. Single-dish monitoring with a dense enough time sampling and a carefully calibrated polarimeter, in combination with 
VLBI observations, are suited to detecting and locating a periodic modulation of the circular maser polarization due to a precess- 
ing Jovian magnetosphere. The phenomenon will be rare, because a favorable arrangement of maser and magnetosphere is needed. 
Otherwise the polarization may be below the detection threshold, especially if the maser is unsaturated. Though exhibiting a qualita- 
tively similar modulation, linear polarization is likely to suffer more from confusion due to dilution of the magnetosphere within the 
maser cross section, even in VLBI observations. 
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1. Introduction 

Most SiO masers are hosted by the extended atmosphere of 
evolved stars of low- to intermediate mass in the upper part of 
the asymptotic giant branch, before they evolve towards central 
stars of planetary nebulae (for a review see e.g., Herwig 2005). 
Under certain conditions, these masers trace the magnetic field 
by virtue of their polarization. After the first theoretical consider- 
ation of a Zeeman laser (Sargent et al. 1967), astronomical maser 
polarization has been extensively discussed, e.g., by Goldreich 
et al. (1977), Western & Watson (1984, for linear polarization), 
Deguchi & Watson (1986, for circular polarization), and Elitzur 
(1991). The models are compared by Gray (2003). The environ- 
ment of SiO masers in AGB atmospheres is characterized by a 
dense wind driven by the stellar pulsations. Its interaction with 
planets engulfed in the stellar atmosphere has akeady been ad- 
dressed by Struck-Marcell (1988) who investigated the hypoth- 
esis that SiO masers may form in the magnetosphere of Jovian 
planets. Struck et al. (2002, 2004) propose scenarios along the 
same lines, and model the dynamic effects of episodic accre- 
tion onto the planets. In a different context, the consequences 
of planets in the stellar atmosphere have been investigated by 
Soker (2001), who shows that Jovian and even Earth-like plan- 
ets can efficiently spin-up the star after migration into the stellar 
atmosphere, such that it later may form an elliptical rather than 
a spherical planetary nebula. 
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Here I propose to resume the idea of Struck-Marcell (1988) 
that SiO masers (though not all) may originate in Jovian mag- 
netospheres, leading to circular maser polarization (Barvainis et 
al. 1987; Herpin et al. 2006), which may be as strong as ~ 10% 
of the Stokes I flux. I show that the densely sampled time se- 
ries of all Stokes parameters yield a rare, but observable polar- 
ization signature characteristic of a precessing magnetosphere. 
Since AGB stars are slow rotators, a precession period as fast as 
~ 10 h hints at a planetary magnetosphere. SiO masers are an 
ideal tool for sounding Jovian planets, since they arise in a zone 
extending from 1.5 to 7 AU distance from the star where they 
are radiatively and collisionally pumped (Gray et al. 2009), be- 
fore the SiO molecule becomes underabundant for maser action 
thanks to condensation in the dust envelope. The SiO maser zone 
thus coiTesponds to a distance from the star where the solar sys- 
tem harbors its giant gas planets (Jupiter is at about 5 AU from 
the sun, and Saturn at about 10 AU). A direct comparison be- 
tween the solar sytem now and in the AGB phase is not valid, 
though, because planetary orbits may be dragged inwards by 
tidal friction, or expelled by the kinematical effect of stellar mass 
loss (Villaver & Livio 2006, further references therein). Arguing 
with Soker (1996) that elliptical planetary nebulae formed from 
stars spun up by inward migration of a planet, it is assumed here 
that the zone in question harbors Jovian planets, since 50 % of 
all planetary nebulae are not elliptical (Soker 2001). 

The first observational evidence of extrasolar planets around 
evolved stars was given by Silvotti et al. (2007) who indirectly 
detected a planet around a post-red giant star, measuring the 
pulse modulations of the variability of the latter As for the plan- 
ets detected around pulsars (with a similar method, measuring 



2 



H.W. Wiesemeyer: The modulation of SiO maser polarization by Jovian planets 



the period variations of the latter) by Wolszczan & Frail (1992), 
it is not yet certain whether these planet-sized bodies survive the 
supernova explosion, or whether they formed in the debris disk 
forming in reponse to a "fallback" (Wang et al. 2006) of ejected 
matter. Most extrasolar planets have so far been detected around 
main-sequence stars by radial velocity measurements (for a re- 
view see Santos 2008). As for the earliest phases of stellar evo- 
lution, a "hot Jupiter" has recently been discovered around the T 
Tauri star TW Hya (Setiawan et al. 2008). Around evolved stars, 
ten companions of substellar mass have been discovered to this 
day, the last one being the KO giant HD 17092 (Niedzielski et 
al. 2007, further references therein). There is no evidence yet 
from AGB stars, mainly because the strong, dense wind makes 
radial velocity measurements of the star difficult. The following 
models of maser polarization modulation due to precessing plan- 
etary magnetospheres are intended as a complement to exoplanet 
searches around AGB stars with different methods. 



2. Description of the model 

2.1. Physical characteristics 

SiO is non-paramagnetic, with a Zeeman splitting of gQ[s"'] ~ 
10^B[G] (cf. Watson & Wyld 2001). For a 1000 K gas and the 
magnetic fields considered here (of up to 10 G), the spectral 
linewidth is greater than the Zeeman splitting by a factor of 
~ 100. This weak Zeeman splitting still produces a fractional cir- 
cular polarization (hereafter pc) of ~ 1 %, which is observable 
for a carefully designed polarimeter (e.g., Thum et al. 2008). 
While linear polarization may not only be caused by magnetic 
fields, but also by anisotropic pumping in conjunction with the 
Hanle effect (Bujarrabal & Rieu 1981; Western & Watson 1984; 
Asensio Ramos et al. 2005), circular polarization is a tracer of 
magnetic flux density as long as the Zeeman splitting is 
much higher than the stimulated emission rate R and the loss 
rate F. As soon as R is close to the circular polarization will 
become intensity-dependent when the maser radiation saturates 
(Nedoluha & Watson 1994), thus ceasing to be a reliable tracer 
of magnetic fields. 

To keep the model simple and to separate the predicted ef- 
fects from those issued by models involving anisotropic pump- 
and loss-rates (see Sect. 3.3), the pump and loss events have been 
assumed to be isotropic here. Following Watson & Wyld (2001), 
the approach of phenomenological maser theory is used. I model 
the polarization properties of the v - 1,7 = 1-0 transition with 
rest frequency = 2.71 x 10" s"', allowing for a straightfor- 
ward analytical formulation of the problem. The rate equations 
yield the population differences «+, «(> and «_ between the split- 
up upper state {J - \,M - -1-1,0,-1, respectively) and the lower 
state (7 = 0, M = 0). These population differences are normal- 
ized by the ratio of the differential pump rate (into the 7=1 
respectively 7 = levels) to the loss rate (assumed to be the 
same for both levels), and are a function of molecular velocity. 
Because of the isotropic pumping by collisions or unpolarized 
radiation assumed here, the pump rates are equal into the sub- 
levels of the 7=1 state. Collisionally induced transitions be- 
tween magnetic sublevels of a given 7 state are neglected. This 
yields 



Table 1. Doppler factors for evaluation of the stimulated emis- 
sion rates 
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Here the stimulated emission rates R± and Rq for the transitions 
with AM = +1 and are normalized by the saturation intensity 
Is - Efia)^r/37Tc^A (A being the Einstein coefficient for spon- 
taneous emission) and are, for a given observing frequency oj, 
evaluated at frequencies a>D, where D are the Doppler factors 
given in Table [T] 

These Doppler factors account for the fact that molecules at 
different line-of-sight velocities - due to their velocity distribu- 
tion and, if present, a velocity gradient (not considered here) - 
couple to different Zeeman components. The stimulated emis- 
sion rates Rj- and Rq are given by (Goldreich et al. 1973) 

R± - I-t(l + cos^ y) + sin^ y(Q-tCos2ri - U±sm2rf) 

+2y±cosy, (4) 
/?o = 2 sin^ y(lQ - Qo cos 2t] + Uo sin 2r]) (5) 

where 7, Q, U and V are the Stokes parameters, normalized by 
the saturation intensity and multipliecO] by An, and the indices 
+, indicate the frequency at which these quantities are to be 
evaluated (Table [T]). Here, y is the angle between the magnetic 
field and the line-of-sight, directed outwards from the drawing 
plane of Fig.[T] rj is the angle between the z axis and the pro- 
jection of the magnetic field onto the plane of the sky (i.e., the 

drawing plane in Fig.[TJ, Bsky = -^By + BI. These angles are 

calculated from 

B^ B^ 
cos y - , cos 7/ = -^^^ . (6) 
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For a static magnetic field, the coordinate system is conveniently 
chosen such that 77 = 0° (i.e., the system in which Stokes U 
vanishes in the absence of magnetorotation). 

2.2. Radiative transfer 

Since the light-travel time is negligible with respect to the time 
scale on which the magnetic field will vary, the stationary radia- 
tive transfer equation is used: 
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' For a linear maser considered here, intensity means its solid-angle 
average multiplied by 4 n. 
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As usual in maser theory, the source term taking spontaneous 
emission into account is insignificant here, and t is the opacity 
along the line of sight for the polarized plus unpolarized radia- 
tion. The Miiller matrix K of the absorption coefficients has the 
elements 



K = 



A B F C\ 

B A E G 

F -E A D 

C -G -D A 
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If the vertical polarization is along the projection of the magnetic 
field onto the plane of the sky, Bsky, the elements of K have the 
following meaning: Coefficients A, B, C and F are the absorp- 
tion coefficients for different polarization states, while D, E and 
G describe anomalous dispersion effects (Landi degl'Innocenti 
& Landi degl'Innocenti 1981). Elements F and G are assumed 
to be negligible, since gO. » R. This condition is necessary here, 
otherwise nothing about the magnetic field can be inferred from 
measurements of pc, since a given molecular state would be de- 
excited by a stimulated emission before a full Larmor precession 
is accomplished. The matrix elements A, B, and C are, for a two- 
level system with split-up upper level, given by (Watson & Wyld 
2001), 
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The terms D and E describe magneto-rotation, i.e., the conver- 
sion of linear to circular polarization, or -in other words - the 
generation of the latter by means other than the Zeeman effect. 
They are given by 
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dr 
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where the extra term Idrj/dr in E accounts for the rotation 
of the coordinate system, while the photon propagates along 
the line-of-sight, such that vertical polarization is always along 
5sky. This extra term is a lengthy expression calculated from the 
magnetic field Eqs. 19 to 21 and the transformation defined by 
Eqs.|22]|23](see below). To correctly handle phase shifts induced 
by magneto-rotation, the profile function / is complex, where 
/l"^' and /I'' are its real and imaginary part, respectively, at the 
normalized frequency offset 
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where Aa>D is the Doppler width of the line, and wo the res- 
onance rest frequency of the transition. Here I use Acj£, - 
1 km s ' , which corresponds to the thermal linewidth of SiO in a 
gas at lOOOK. 

Correspondingly, and are the real and imaginary parts 
of the profile function at 

y- = (w- [wo -gQ/2]) 
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and f^'^ and /® at 
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The complex profile function f(a>) is given by 
<i>{v)dv 
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where Y is the decay rate of the excited state (5 s"' for the SiO 
masers considered here, see Nedoluha & Watson 1994), and <l)(v) 
is the velocity distribution of the masing gas, which is taken to be 
Maxwellian. The velocity distribution of molecules in different 
substates differs from a pure Maxwellian distribution, because 
Eqs. 1-3 are velocity-dependent. The real and imaginary parts of 
fico) are known as the Voigt function 'H and the Faraday- Voigt 
function T, respectively, namely 



f(v, fl) = ^ YH{a, v) + 2/ T(a, v)] 



(18) 



with a - r/AwD. For '7/ and T I used the rational approximation 
given by Humlicek (1982), as provided in an optimized form by 
Schreier(1992). 

Without a detailed model at hand, I include neither the mag- 
netopause nor the magnetotail, but assume that the SiO maser 
slab crosses a magnetic dipole field (in spherical coordinates, de- 
fined by the symmetry axis of the dipole field, e.g., Lewis 2004), 
given by 
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where M - 34Grj is the magnetic dipole moment (rj is Jupiter's 
mean radius). The justification of using such a magnetic field 
component is discussed in Appendix A. We now transform the 
magnetic dipole field from its stationary frame into a coordinate 
system defined by the planetary rotation axis, Cz = (0,0, 1)^, 
the line-of-sight towards the observer, Cx - (1,0,0)^ and Cy = 
(0, 1 , 0)^, such that (Cx, Cy, Cz) form a left-handed coordinate sys- 
tem (Fig.[TJ. The transformation is defined by the Euler angles 
O - Int/T = D.t (where T is the rotation period, and the incli- 
nation of the magnetic dipole axis against the rotation axis), and 
is obtained by 



(22) 
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with 
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sin sin Qx - sin cos Qx cos , 



(23) 



(in Fig.[ri O = 90° has been chosen). For T and 0, values close 
to those for Jupiter are used, i.e., T = lOh and - 10°, as- 
suming that the planet is still too far away from the star to lose 
angular momentum to spin up the latter (Soker 2001). Keeping 
the typical values for giant gas planets in the solar system is 
therefore not too farfetched in this context. 



3. Results and discussion 

In the following, the model results (Figs.|2]to |7]i are shown as 
grey-scale plots (fractional polarization or polarization angle vs. 
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Fig. 1. Sketch of the model with magnetic field lines. The x,y,z coordinate system (solid lines) is defined by the rotation axis of 
the planet (z axis) and the line of sight to the observer (x axis, directed from the drawing plane outwards). The coordinate system 
(dashed lines) of the magnetosphere is defined by the tilt angle of the dipole axis with respect to the rotation axis. The 
rotation phase here is chosen such that the x' axis is identical to the x axis. The cross section of the planet is depicted by the black 
circle, the maser slab by the grey one (for models lb, 2b). The location of the Chapman-Ferraro magnetopause is indicated by the 
dotted line. The top right insert shows the (x, y) plane with the location of the magnetopause (dashed line) and the maser slab as a 
grey rectangle (field lines omitted). The bottom right insert shows a zoom onto a particular line of sight, marked by a cross, within 
the SiO maser bundle. The arrow indicates the local magnetic field direction and the definition of angle r\. 

Table 2. Summary of models 
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phase and velocity) and as polarization spectra at a given ref- 
erence phase of the planetary rotation, or as a time series at a 
given reference velocity. To search for strong magnetic fields, 
maser features with a strong fractional circular polarization (pc) 
are needed (which also suff'er less than linear polarization from 
the confusion in the main beam of the spatially unresolved ob- 
servations). I therefore use the phase and velocity of the maser 
feature with the strongest po, as reference. 



3A. Saturated maser 

The first model la (see Table 2 and Fig.|2]i is for a maser with 
a pencil beam in the planet's equatorial plane, at 3 rj from the 
center of the planet, with a length of 10 rj. The opacity is scaled 
such that the peak intensity is 10 /s, i.e., the maser is saturated. 
Because of both the small tilt of the magnetic dipole axis against 
the rotation axis, and the absence of a toroidal magnetic field 
component B^, cos y ^ 1 and the linear polarization does only 
slightly exceed the canonical value of 33 % (Goldreich et al. 
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Fig. 2. Saturated SiO maser (/peak = 10 /s), pencil beam at 
r - 3rj, in the planet's equatorial plane (model la). Inserts: spec- 
tra and time series at the phase (indicated by a horizontal dashed 
line), respectively velocity (indicated by a vertical dashed line), 
of maximum circular polarization. Each dot represents a point 
of the time sampling. Top: fractional linear polarization (p^ vs. 
velocity and phase (contour levels 15 % to 33 % by 3 %). Center: 
position angle x of linear polarization East from North, contour 
levels -18° to +18° by 4°). Bottom: fractional circular polariza- 
tion (pc, contour levels from -1.2 % to +1, 2 % by 0.2 %, zero 
contour suppressed). 



1973) for a saturated maser propagating perpendicular to the 
magnetic field. (Our excess stems from line-of-sight elements 
with cos^ -y > 2/3.) The peak-to-peak variation in the fractional 
linear polarization, pL, is 6.6% at line center, and 3.4% at the 
velocity of maximum circular polarization. The circular polar- 
ization displays the archetypical S-shaped line profile of an un- 
resolved Zeeman pattern whose peak-to-peak amplitude varies 
in time by 2.6 %. This effect will be observable provided that the 
maser is sufficiently strong. The relatively strong circular polar- 
ization is not only produced by the Zeeman effect (term C in 
Eq.[8]), but also increased by the rotation of the magnetic field in 
the maser amplification zone along the line-of-sight, rotating the 
linear polarization from Stokes Q to Stokes U via the E-term in 
Eq.m and then its conversion from Stokes U to Stokes V via the 
D-term. 

In observational practice, the situation is often less clear, for 
two reasons. First, while polarization monitoring with single- 
dish telescopes offers the advantage of being free of fluctuations 
due to baseline synthesis effects in VLBI observations, it lacks 
the spatial resolution of the latter What is obtained is instead a 
spectral blend of Zeeman features from individual maser spots, 
centered at different velocities. If the full width to zero Stokes 
V from an individual feature is below the velocity separation of 
two spatially separated maser features, observations with a high 
enough spectral resolution will help. If not, this results in confu- 
sion. Even the spectra of spatially resolved maser spots may be 
concerned in the case of a velocity gradient along a single maser 
slab, leading to velocity redistribution. As akeady stated above, 
the latter case is not handled here, with the rationale of restrict- 
ing the models to the basic cases before dealing with the com- 
plications (which, in the limit of a uni-directional slab maser, 
can of course be handled by an extension of the theory described 
above). Nature may help us here, because the magnetic flux den- 
sity in a planetary magnetosphere may be enhanced with respect 
to that of a magnetic field of stellar origin, but this hope needs to 
be ascertained. We cannot expect that all detections of a strong 
Stokes V in SiO masers (corresponding to a few up to 20 G, 
Herpin et al. 2006) are owing to planetary magnetospheres. 

The second caveat concerns both VLBI and single-dish ob- 
servations. In reality, SiO masers in Mira stars have a typical size 
of at least ~ 10'^ cm (constrained both from VLBI observations, 
e.g.. Philips & Boboltz 2000, and models, BujaiTabal 1994), 
which is much larger than the length scale on which the magnetic 
dipole field varies in model la. To simulate the effects of dilut- 
ing the magnetosphere in the maser cross section, a cylindrical 
maser slab (model variant lb) has also been modeled as a collec- 
tion of maser rays with a bundle diameter of 3 o (i.e., two orders 
of magnitude below observed spot sizes). This corresponds, at a 
distance of e.g., 230 pc, to an angular size of 0.01 mas (i.e., still 
unresolved by VLBI observations). Assuming a larger and more 
realistic size would require including the magnetopause and the 
magnetotail in the model, and possibly the detailed dynamics 
of mass accretion onto the planet, which clearly is beyond the 
scope of this paper. I therefore assume that the polarization sig- 
nals modeled here are on top of a linear or circular polarization 
background, which does not vary on a timescale comparable to 
that of the planetary rotation. We know that the lifetime of SiO 
masers is unlikely to be longer than a cycle of the stellar vari- 
ability, because the features are disrupted by shock waves driven 
by the stellar pulsation (Humphreys et al. 2002). Multiepoch 
VLBA imaging of SiO v = 1,7 = 1-0 masers impressively 
confirmed this picture (the "movie" of TXCam, Diamond & 
Kemball, 2003). Long-term polarization fluctuations may occur 
on shorter timescales (several weeks at least, Glenn et al. 2003). 
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They come from slow readjustments of the magnetic field. These 
timescales are long enough to monitor at least 10 periods of the 
planetary rotation, allowing us to separate the short-term polar- 
ization fluctuations modeled here from the long-term ones of po- 
larized and total maser flux. 

The polarization background is not necessarily dominant, 
even in single-dish observations. Linear polarization is reduced 
by the combined action of both the predominantly tangential po- 
larization of the maser spots distributed along a ring (providing 
indirect evidence of radiative pumping, Desmurs et al. 1999) and 
cancellation of polarization vectors in the observing beam. As 
for circular polarization, we may have reasons to expect that it is 
enhanced towards the lines of sight crossing the strongest mag- 
netic fields, i.e., close to the planet, as suggested above. Any ob- 
servational attempt to discover the features modeled here will 
therefore focus on measurements of the circular polarization. 
Fig.|3]shows that the maximum linear polarization is reduced by 
about 5.3 % due to dilution (even in VLBI observations, since 
they do not resolve the maser spot either), and incoherent mix- 
ing correspondingly reduces the peak-to-peak variation of the 
polarization angle. As expected, the circular polarization basi- 
cally remains the same. (Because of the mixing of maser rays 
with smaller Zeeman splittings, the velocity of maximum pc is 
reduced.) At the velocities of maximum circular polarization, the 
relative peak-to-peak variation of the linear polarization has dou- 
bled with respect to the pencil-beam model (la), thanks to the 
shift in the velocity of maximum circular polarization (which is 
used as reference). 

3.2. Unsaturated maser 

The saturation of masers is diflicult to determine experimentally, 
mainly due to its dependence on the observationally poorly con- 
strained beaming angle (Watson & Wyld 2001). The polarization 
of unsaturated masers is weaker than that of their saturated coun- 
terparts, because only the latter convert each pumping event into 
a maser photon. A fully unsaturated maser (F <K R) is unpolar- 
ized if both the seed and pump radiation are unpolarized. The 
model 2a (Table 2 and Fig.lUi is for a pencil beam with 0.1 /s, 
and yields a maximum linear polarization of 4.4 % at line center 
(with a peak-to-peak variation of 1.3 %), which drops to 2.6 % 
at the velocity of strongest circular polarization (with a peak-to- 
peak amplitude of 0.8 %). As expected in unsaturated masers, 
the line profile of Stokes V is identical to that of thermal radia- 
tion (Fiebig & Gusten 1989, see also Watson & Wyld 2001), and 
the linewidths are narrowed with respect to the thermal width 
(cf. Watson & Wyld, Fig. 2). Therefore, the circular polarization 
is drastically reduced (to pc = ±0.009 %) for a maser ray in the 
equatorial plane, and is no longer detectable with current sen- 
sitivity. Again, a cylindrical maser slab (model 2b, Fig.|5] elec- 
tronic version) does not drastically alter the result for circular 
polarization, but reduces the linear polarization at the reference 
channel to a maximum of 1.5 %, with a peak-to-peak variation 
of 0.3 % (i.e., mean value and amplitude are scaled down by a 
factor 1 : 2). 

For the detectability of the circular polarization of an unsat- 
urated SiO maser from a Jovian magnetosphere, its location in 
the latter is important. Because the circular polarization here is 
proportional to the line-of-sight component of the magnetic field 
(hereafter Bios), it is naturally enhanced for maser rays crossing 
the rotation axis, which have a maximum Bios before and after 
the crossing. This is demonstrated by model (2c) in Fig.|6](elec- 
tronic version), where the maser ray is otherwise at the same dis- 
tance from the planet as in model (2a). The peak-to-peak vari- 
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Fig. 3. Same as Fig.|2] but for a bundle of SiO masers, with a 
slab diameter of 3 o (model lb). For easy comparison, the plot 
limits for the inserts are the same as those for the equivalent 
pencil beam model (Fig.|2]i. Top: pi^, contour levels from 15 % 
to 27 % by 3 %. Center: x, contour levels from -10° to -Hi 0° by 
5°. Bottom: pc, contour levels from -1.2 % to H- 1.2 % by 0.2 % 
(zero contour omitted). 



ation of pc — ±0.29 % now becomes detectable with a care- 
fully tuned polarimeter (whereas the maximum linear polariza- 
tion drops to PL = 1.5 %). 
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For lines of sight close to the rotation axis, B^ty changes 
strongly as a function of both velocity and rotation phase. Model 
2d (Fig.|7] available online) is for a pencil beam 3 rj above the 
equatorial plane, and slightly offset (0.5 rj) from the rotation 
axis. The dependency of the polarization angle on both Bios and 
the magnetic field component projected onto the plane of the sky, 
Bsky leads to a complicated pattern. The polarization angle flips 
by 90° are because of the critical cos^y = 1/3, which decides 
whether the linear polarization is along B^ky or perpendicular 
to it (Goldreich et al. 1973), and of the variable maser satura- 
tion across the line profile (hence the strong spectral variations). 
Observationally, these flips have been proven by VLBI polariza- 
tion maps (i.e., Kemball & Diamond 1997) and hint at variations 
of the magnetic field either due to a planetary magnetic field as 
modeled here or a magnetic field component of stellar origin. 

3.3. The impact of refinements of fhe radiative model 

3.3.1. Anisotropic pumping 

SiO masers are believed to be radiatively pumped by infrared ex- 
citation of vibrationally excited states (Gray et al. 2009, further 
references therein). The anisotropic distribution of the pump ra- 
diation leads to an unequal population of the magnetic sublevels 
and therefore enhances the linear polarization (Elitzur 1996), 
which generates a circular polarization in excess of its value for 
isotropic pump rates if magnetorotation is present (as in the mod- 
els presented here). This bkefringent conversion from linear to 
circular polarization would make a detection of fluctuations of 
PC possible even in the unsaturated regime. The same presum- 
ably holds for anisotropic loss rates. In the case of anisotropic 
pumping, the critical value of y for producing linear polarization 
parallel or perpendicular to Bsky will not be the 55° of Goldreich 
et al., though the linear polarization will still have one of these 
orientations with respect to the magnetic field. 

3.3.2. More complex maser geometries 

A refined model would imply not only a more adequate descrip- 
tion of the magnetosphere crossed by the SiO maser, but also 
a maser geometry more complicated than the uni-directional 
cylindrical maser slabs considered here. The first step to a 
more realistic description of the maser action is to consider bi- 
directional masers (i.e., two opposed propagation directions in 
a linear maser). For a constant magnetic field, Watson & Wyld 
(2001) show that the circular polarization remains nearly unal- 
tered between uni- and bi-directional masers, whereas the linear 
polarization has a tendency to be increased in the last (while 
qualitatively keeping its dependency on maser saturation and 
cos y). 

3.3.3. Higher rotational quantum numbers 

The comparison between the v = 1,7 = 2-1 maser transi- 
tion and itsv - I, J = 1-0 counterpart (both taken to be 
uni-directional) shows that their circular polarization is not dras- 
tically different. As for linear polarization, Deguchi & Watson 
(1990) found that the polarization of the 7 = 2-1 transition 
is actually suppressed with respect to the 1-0 transition when 
(gQ.)^/r » 7J » gQ., while observationally the former tends to 
be enhanced (Mcintosh & Predmore 1993), such that strong lin- 
ear polarizations can only be achieved by other than magnetic 
means (e.g., the aforementioned anisotropic pumping, as pro- 
posed, in this context, by Watson & Wyld 2001). However, this 




Fig. 4. Unsaturated SiO maser (/peak = 0.1 /s), pencil beam 
at r = 3 o, in the planet's equatorial plane (model 2a). The 
linewidths are narrowed with respect to the thermal width, as 
expected for unsaturated maser action. Top: pi^, contour levels 
from 0.6 % to 4.2 % by 0.6 %. Center: x, contour levels -12° 
to 12° by 3°. Bottom: pc, contour levels from -0.009% to 
-hO.009 % by 0.002 %. 



regime does not allow us to analyze the magnetic field geom- 
etry by means of polarization measurements, because a single 
Larmor precession of the molecule is destroyed by a stimulated 
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emission before its completion. I therefore did not consider this 
case here. 

3.4. The impact of instrumental limitations 

3.4.1. Spatial resolution 

While the quasi-stationary linear polarization tends to decrease 
due to mutual cancellation within the beam of spatially unre- 
solved observations, the variations in the circular polarization 
from different maser spots would need to fluctuate in anticor- 
relation in order to cancel. This is highly unhkely, even if all 
variations were caused by the effect discussed here, and even 
more if other phenomena - e.g., jets or coronal flux loops, oc- 
curring on different timescales - were involved. These scenarios 
are discussed in Wiesemeyer et al. (2009). The observational re- 
sults of that work, and subsequent polarization monitoring from 
December 2008, indicate that fluctuations of the Unear and cir- 
cular polarizations may appear simultaneously (as shown here) 
or separately (in contrast to the models presented in this paper). 
The possibility of a blend of several unrelated phenomena calls 
for a coordinated observational effort as proposed in our 2008 
paper. The direct imaging of planetary wake flows in the atmo- 
spheres of the nearest AGB stars with ALMA will definitely be 
an asset. According to Cherchneff (2006), we can expect suitable 
molecular high-density tracers to be sufficiently abundant in the 
atmosphere, thanks to the non-equihbrium chemistry dominated 
by shock propagation. 

3.4.2. Spectral resolution 

The Stokes spectra shown in Figs. 2 to 7 are calculated with 
a frequency resolution of 27rAv - 10 kHz, corresponding to a 
velocity resolution of 0.01 kms~\ which in practice cannot be 
achieved with existing correlation polarimeters. After smooth- 
ing with a box-like kernel, 12 frequency steps wide, the circular 
peak polarization of model la (Fig. 2) drops from ±1.27% to 
±0.94%, which is still observable (i.e., the spectral smoothing 
does not lead to an order-of-magnitude loss of circular polar- 
ization). The decrease in circular polarization due to the Umited 
spectral resolution depends on the width of the Zeeman feature. 
A comparison between the pc spectra in Figs. 2 to 7 shows that 
the relative decrease of the circularly polarized signal remains 
comparable for the unsaturated masers. However, here the weak- 
ness of the polarization will make a detection difficult, if not im- 
possible. The polarization monitoring therefore has to be done 
with the best possible spectral resolution. 

4. Conclusions 

I have shown that densely sampled polarization monitoring is 
able to reveal precessing Jovian magnetospheres engulfed in the 
atmospheres of AGB stars by a periodic modulation of the po- 
larization characteristics of SiO masers crossing them. The ex- 
pected period is that of the rotation of giant gas planets (i.e., 
~ 10 h if the rotation period did not change with respect to typ- 
ical solar system values), which allows us to distinguish such a 
fingerprint from long-term variations of the polarized (Glenn et 
al. 2003) and unpolarized maser flux (Humphreys et al., 2002; 
Diamond & Kemball, 2003), which occur on timescales at least 
an order of magnitude longer than those considered here. As for 
the distinction from pseudo-periodic polarization variations due 
to stellar magnetic flux loops or magnetic clouds, everything de- 
pends on the speed of the SiO masers relative to the latter. At a 



speed of 1 km s" , a maser covers a distance of 5 rj within 1 h, 
i.e., the patterns modeled here are only visible if the maser speed 
across the planetary magnetosphere is < 10 km s~\ or if the pre- 
cession period is < 10 h. Otherwise, polarization fluctuations are 
expected that not only contain temporal variations, but also spa- 
tial ones, and the patterns will be either compressed or stretched. 
A planetary magnetosphere exposed to an AGB wind has its 
magnetopause closer to the planet than e.g., Jupiter's magne- 
topause, which is exposed to the faster, but less dense solar wind. 
As a matter of fact, the case of a planetary magnetosphere in 
an AGB wind may be closer to that of the terrestial magneto- 
sphere, whose magnetopause is as close as a few Earth radii, and 
can even be temporarily disrupted in case of strong solar storms. 
Likewise, the magnetic dipole component of a planetary mag- 
netic field engulfed by the atmosphere of an AGB star can be 
expected to be quite compact (not more than 10 planetary radii). 
To observe the modeled features, an SiO maser has to cross it, 
and either needs to be saturated to produce a measurable circu- 
lar polarization, be unsaturated if the pump and loss rates are 
sufficiently anisotropic, or has to cross a part of the magneto- 
sphere with a strong line-of-sight component of the magnetic 
field. This leads to the conclusion that the phenomenon will be 
rare and, if detected, only be seen in a narrow range of velocities, 
and preferentially in circular polarization-enhancing regions of 
strong magnetic flux density (whereas linear polarization tends 
to cancel out if the polarization angle strongly varies across the 
observing beam). Furthermore, the dynamics of the gas close to 
the planet possibly accreting matter needs to be such that the 
line-of-sight velocity coherence allows for the maser action. The 
abundance of SiO molecules required to form a strong enough 
maser can be expected to be less of a problem. The evaporation 
of Galilean moons has been proposed as an additional reservoir 
for the gas-phase SiO (Struck-Marcell, 1988). 

Recently, Wiesemeyer et al. (2009) have detected a pseudo- 
periodic variation of the fractional circular polarization towards 
two Mira stars, R Leo, and V Cam, which may hint at precess- 
ing magnetospheres in their atmospheres. Several periods will 
have to be monitored, to better distinguish the phenomenon from 
other rapid, but transient magnetospheric events. The case of 
R Leo is especially intriguing, since VLBA maps of the 43 GHz 
SiO maser (Cotton et al. 2004, 2008; Soria-Ruiz et al. 2007) 
show elongated SiO features pointing radially away from the 
star, reminiscent of jet-Uke features. The velocity information 
from the observed polarization fluctuations, together with the 
astrometry of such elongated features, suggest orbital motion in 
at least one of them, and strengthens the case for a planetary 
wake flow as simulated by Struck et al. (2004). Dense polariza- 
tion monitoring of these features may be especially rewarding, 
and a spatially fully resolved polarization monitoring will cer- 
tainly be highly conclusive. However, VLBl techniques will suf- 
fer from a confusion between polarized flux variations and vari- 
ations due to Earth rotation synthesis of visibility components 
(the timescale for both is unfortunately the same). This calls for 
a combination of single-dish polarization monitoring with VLBl 
imaging of the Stokes 1 emission. The cross-identification of po- 
larization fluctuations with VLBl features is possible if there 
is, in the single-dish observations, no blend of spatially sepa- 
rated SiO maser spots at the same velocity or if the VLBl fea- 
tures show an unambiguous association with a planetary wake 
flow. Detection of orbital motion in such features will provide 
a touchstone for the scenarios presented here. The aim of this 
paper has been to motivate coordinated observational efforts, so 
far unprecedented, and to show that meaningful results can be 
expected to emerge from them. If so, this will call for advanced 
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modeling. The work published here should therefore be consid- 
ered as the initial step in an iterative refinement. 
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Appendix A: The magnetosphere model 

Before modeling a Jovian magnetosphere exposed to an AGB 
wind, it has to be shown whether it can resist the latter, i.e., 
whether it is able to form a magnetopause far enough from the 
planet to allow for a protected magnetic dipole field, such that 
the ram pressure of the inflowing wind equals the magnetic pres- 
sure of the magnetosphere (Chapman-Ferraro boundary, see e.g.. 
Kennel, 1995). This is the case if the magnetic flux density in the 



magnetopause, B^^, is 



(lO-K'gcm-^) (lOkms-i)^' ^^'^^ 

where pw and vw are the density and velocity of the wind, re- 
spectively (Struck-Marcell, 1988), scaled to the typical values 
for an AGB wind. The assumed Jovian magnetic dipole field 
reaches this value in the equatorial plane at a distance 



B 



0.05 



M 



(A.2) 



For a magnetic moment of M[G] - 34 (eight times the Jovian 
value), a dipole magnetosphere can extend to up to 8.8 Jovian 
radii (substellar distance from the planet). For a weaker mag- 
netic field, the magnetopause would be closer to the planet, but 
it would still be possible to have an SiO maser crossing the pre- 
cessing magnetosphere in an extended magnetotail. 

A comparison with the solar system may be instructive here. 
The Jovian magnetopause extends to 92 rj, which is about dou- 
ble what would be expected from Eq. lA.ll for solar system pa- 
rameters (pw = O.lAcmT^, v'w = 400 kms"'), showing that the 
Jovian magnetosphere is more complex than a magnetic dipole 
field protected by a magnetopause (Alexeev & Belenkaya 2005). 
Only a strong enough interplanetary magnetic field can push the 
magnetopause beyond this limit. In the atmosphere of an AGB 
star, even more complex phenomena may be expected, because 
of the high density of the wind and the possibility of a reconnec- 
tion between the planetary and stellar magnetic field. For com- 
parison, the terrestial subsolar magnetopause may be as close at 
6 r-E to the Earth's center, depending on space weather (Shue et 
al. 1997), i.e., at a distance (relative to the Earth's radius) rather 
similar to that of the Jovian model planet (relative to rj) than 
to that of Jupiter. Likewise, the location of a planetary magne- 
topause in an AGB wind depends on the variations in the latter, 
especially in response to the stellar pulsations. These timescales 
are, however, far longer than those considered here (~ lOh). To 
avoid a perturbation of a precessing magnetosphere by the AGB 
wind and the accretion of matter onto the planet, the magnetic 
flux needs to be decoupled from the matter via ambipolar diffu- 
sion. As for the latter, I estimate its speed Vad following Hartquist 
& Dyson (1997), but using values closer to the magnetosphere 
modeled here, by 



Vad 



lOOOkms-' i-^ii-^Vi^] 
llG/ \0.1rj/ \\Q-(>I 

\10'«cm-3j 



(A.3) 



where L ^ 0.1 rj is the typical length scale on which the mod- 
eled magnetic field varies, B ^ 1 G is the typical magnetic flux 
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density in the model magnetosphere, Xi ^ 10"^ the ionization 
fraction, and nu ^ lO'^gcm"-^ the expected hydrogen density 
in the maser slab. The major uncertainty in Eq. lA.3l comes from 
the gas density (especially if shocks are involved) and the ion- 
ization fraction. In any realistic case, the suggested ambipolar 
diffusion speed of ~ 1000 km s"' is much faster than the pre- 
cession speed of the magnetosphere (2nrs/T ~ 12kms"' with a 
period of T - 10 h), and the magnetic flux and matter can decou- 
ple efficiently, such that the magnetic field lines diffuse through 
the maser slab. This conclusion holds for the localized magnetic 
field considered here. As for the global magnetic field pervad- 
ing the AGB atmosphere, Soker (2006) rules out any dynamical 
importance, because of the weak coupling between the magnetic 
flux and the matter, based on a similar argument. 1 conclude that 
a planetary magnetic dipole field can persist in the extended at- 
mosphere of an AGB star. Ohmic dissipation of the magnetic 
flux can be neglected here, otherwise the substantial circular po- 
larizations of SiO masers could not have been observed. 



Appendix B: Numerical calculations 

Because the growth rates of the Stokes parameters in the prop- 
agation of maser photons differ by orders of magnitude at dif- 
ferent velocities, Eq. 7 is a stiff set of differential equations, and 
numerical solutions using Runge-Kutta or multiple-order meth- 
ods are inefficient, even with adaptive step sizes. Because the 
Stokes parameters at different frequencies are not independent 
of each other due to the Zeeman splitting, the calculations have 
to be done along the full zero-to-zero width of the spectral line, 
with a sampling dense enough to resolve the Zeeman feature. 
For a given time step t, the calculations proceed along the line- 
of-sight in steps that are short enough to solve the equation of 
radiative transfer using the matrix exponential function exp (Kr) 
with piecewise constant coefficients, i.e., for a step from grid 
point j to 7 + 1 

Ij+i =exp(KT)Ij (B.l) 

where the elements of the vector Ij are the Stokes parameters en- 
tering the grid cell. Spontaneous emission is neglected, as usual 
in maser theory. Landi degl'lnnocenti & Landi degl'lnnocenti 
(1981) have provided an elegant analytical solution for the ma- 
trix exponential function in closed form. At each sampling point 
along the line-of-sight, we calculate the magnetic field from 
Eqs. 19 - 23, the corresponding angles y and rj from Eq.|6l the 
resulting Zeeman splitting gQ., and the corresponding matrix el- 
ements A to E. For a given line-of-sight element and frequency, 
the stimulated emission rates in Eqs. 1 to 5 are determined from 
the normalized Stokes intensities entering it by a Spline interpo- 
lation of the spectral line profiles in (/, Q, U, V)^. 1 assume a lin- 
ear, unidirectional maser, with unpolarized continuum seed radi- 
ation. For each rotation phase considered, results are shown for 
both a maser spot much smaller than the length scale on which 
the magnetic field typically varies, and for a maser spot consist- 
ing of a number of independent, parallel rays (hereafter called 
bundle). As already mentioned, for distances farther than 8 rj 
away from the planet the model would not be valid anymore, 
and would have to be modified by the shocked AGB wind, the 
magnetopause, and the magnetotail. Without a detailed magneto- 
hydrodynamical model at hand, we assume that all masers of the 
bundle are centered on an impact parameter of 3 o (such that 
the maser is confined within the magnetic dipole field), have the 
same line-of-sight velocity, and are characterized by a total gain 
length of lOrj. The tidal interaction of the masing gas with the 



possibly matter accreting planet leads to a loss of velocity coher- 
ence and thus of maser saturation. I therefore model the effects 
of an incoherent velocity field by comparing unsaturated maser 
action with a saturated one. For the case of a homogeneous mag- 
netic field the code has been successfully tested by comparing its 
results whith those of Watson & Wyld (2001). 
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Fig. 5. Same as Fig. H] but for a bundle of SiO masers, with slab 
diameter of 3 rj (model 2b). Top: pL, contour levels 0.3 % to 
2.7% by 0.3%. Center: contour levels from -15° to +12° 
by 3°. The "noise" in the spectrum of is a numerical arti- 
fact. Bottom: pc, contour levels from -0.007 % to +0.007 % by 
0.002%. 



Fig. 6. Same as Fig.|4] but 3 o above the planet's equatorial plane 
and towards its rotation axis (model 2c). Top: pi^, contour levels 
0.3 % to 1.5 % by 0.3 %. Center: x, contour levels from -16° 
to +160° by 16°. Bottom: pc, contour levels from -0.28 % to 
+0.28 % by 0.04 %. 
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Fig. 7. Same as Fig.|6] but slightly offset from the rotation axis 
(by 0.5 rj, model 2d). Top: pL, contour levels 0.1 % to 0.6 % by 
0.1 % (the small jumps in the insert showing the time series of pL 
are a numerical artifact). Center: contour levels from -16° to 
+ 160° by 16°. Bottom: pc, contour levels from 0.02 % to +0.1 % 
by 0.02 %. 



